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Abstract 


The  concept  of  the  wave  packet  is  used  to  obtain  some  general  results  on  the 
propagation  of  signals  in  dispersive  media.  The  dispersion  of  a  Gaussian  carrier 
pulse  and  a  square  wave  carrier  pulse  in  an  isotropic  plasma  is  carried  out  using 
the  wave  packet  concept.  A  general  analysis  of  transient  wave  propagation  in 
isotropic  plasmas  is  given  using  Laplace  transform  methods.  Solutions  are  given 
in  terms  of  series  solutions  which  may  be  expressed  in  terms  of  Lommel  functions 
Integral  solutions  which  can  be  easily  evaluated  numerically  are  derived  using  the 
convolution  theorem  and  using  contour  integration  techniques.  The  solutions  arc 
useful  for  short  dispersion  lengths. 
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Transient  Signal  Propagation  in 
Lossless,  Isotropic  Plasmas 

Volume  I 


1.  INTRODUCTION 


The  study  of  propagation  of  transient  signals  in  dispersive  media  dates  back  to 
the  early  years  of  this  century.  Following  Einstein's  publication  of  his  special 
theory  of  relativity,  concern  arose  over  the  fact  that  in  regions  of  anomalous 
dispersion  the  group  velocity  is  greater  than  the  free-space  velocity  of  light,  c. 
Since  it  was  believed  that  the  group  velocity  was  the  velocity  at  which  energy  is 
transf  or  ted  by  the  wave,  this  condition  of  anomalous  dispersion  appeared  to  violate 
Einstein's  theory. 

This  paradox  was  correctly  explained  by  Sommerfeld  (1914)  who  showed  by 
making  a  high  frequency  expansion  that  the  very  first  part  of  a  signal,  called  the 
signal  wavefront,  arrives  at  a  given  point  with  the  velocity  of  !ight,  c.  Sommer- 
feld'a  solution  is  valid  for  only  a  short  time  after  the  arrival  of  the  signal  wave- 
front.  By  using  a  saddlepoint  method  of  integration  Brillouin  (1914)  found 
solutions  which  are  valid  in  a  certain  time  interval  following  the  Sommerfeld 
region.  The  signal  in  this  region  is  called  a  precursor  since  it  precedes  the 
arrival  of  the  main  signal.  This  work  of  Sommerfeld  and  Brillouin  has  been 
summarized  in  a  book  by  Brillouin  (1960)  in  which  some  of  the  important  early 
papers  have  been  reprinted. 
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The  arrival  of  the  main  signal  follows  the  precursors.  However,  the  standard 
saddlepoint  method  of  integration  cannot  be  used  in  this  region  since  the  saddle- 
point  is  approaching  a  pole  in  the  complex  plane.  Methods  for  appropriately 
modifying  the  saddlepoint  method  under  these  conditions  have  been  discussed  by 
Cerrillo  (1950),  van  der  Wearden  (1950),  *uid  Clemmow  (1950).  These  techniques 
have  been  used  by  Pearson  (1953)  to  describe  transient  propagation  in  acoustic 
waveguides,  which  have  the  same  dispersion  equation  as  a  lossless  isotropic 
plasma.  Electromagnetic  waveguides  also  have  the  same  dispersion  equation  and 
the  saddlepoint  method  of  integration  has  been  applied  to  these  problems  by 
Cerrillo  (1948),  Namiki  and  Horiuchi  (1952),  and  Karbowiak  (1957). 

The  solutions  discussed  so  far  are  approximate  solutions  that  are  valid  only 
after  the  signal  has  propagated  a  long  distance  through  the  dispersive  "'ed;a.  The 
detailed  consideration  of  these  types  of  solution  will  form  the  subject  matter  of  a 
second  and  subsequent  report  on  this  topic. 

The  present  report  consists  of  two  main  parts.  In  Section  2  the  well  known 
concepts  of  a  wave  packet  and  group  velocity  are  applied  to  the  propagation  of 
pulses  in  a  dispersive  medium.  This  method  has  been  used  to  describe  the  propaga¬ 
tion  of  a  rectangular  pulse  in  a  waveguide  by  Cohn  (1952),  Elliott  (1957),  Knop  and 
Cohn  (1963),  and  Wanselow  (1962).  The  propagation  of  a  Gaussian  pulse  in  a  wave¬ 
guide  using  this  method  has  been  described  by  Forrer  (1958).  These  solutions  only 
apply  to  quasi-monochromatic  signals.  In  Section  3  more  general  solutions  which 
are  valid  for  arbitrary  signals  which  propagate  only  a  short  distance  through  the 
plasma  will  be  considered. 

Jaeger  and  Westfold  (1949)  have  considered  both  the  propagation  of  an  initial 
disturbance  in  a  plasma  when  the  spatial  dependence  is  prescribed  as  well  as  the 
disturbance  that  is  radiated  when  the  time  variation  of  the  electric  field  is  pre¬ 
scribed  at  some  position.  They  have  written  solutions  in  terms  of  convolution 
integrals  and  have  studied  the  Fourier  spectra  of  several  initial  disturbances.  They 
have  applied  their  results  to  the  solar  corona  in  an  attempt  to  explain  some  of  the 
phenomena  associated  with  bursts  of  solar  noise.  Some  numerical  results  based 
on  convolution  integrals  will  be  discussed  in  Section  3.1.  Some  of  these  results 
have  been  given  by  Case  (1965). 

Exact  integral  solutions  may  be  obtained  by  contour  integration.  Cerrillo  (1948) 
has  discussed  such  solutions,  and  in  Section  3.3  two  possibilities  will  be  described 
in  detail  and  numerical  results  will  be  presented. 

For  the  case  of  a  turn-on  sine  wave  the  exact  solution  may  be  written  as  a 
series  of  Bessel  functions.  These  solutions  have  been  discussed  by  Cerrillo  (1948), 
Rubinowicz  (1950),  Gajewski  (1955),  Kovtun  (1958),  and  Knop  (1964).  Cerrillo  (1948) 
and  Kovtun  (1958)  have  expressed  these  solutions  in  terms  of  Lommel  functions. 
These  series  solutions  will  be  discussed  in  Section  3.2. 
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Finally,  the  experimental  work  in  which  observations  of  transient  dispersive 
effects  have  been  observed  should  be  mentioned.  As  has  already  been  pointed  out, 
the  dispersion  equation  is  the  same  for  acoustic  waves  in  fluid-filled  tubes, 
electromagnetic  waves  in  waveguides,  and  electromagnetic  waves  in  lossless, 
isotropic  plasmas.  Each  of  these  situations  has  been  investigated  experimentally. 

The  transient  response  of  sound  pulses  propagating  through  fluid-filled  tubes 
has  been  observed  by  Anderson  and  Barnes  (1953)  and  by  Proud,  Taniarkin,  and 
Kornhauser  (1956).  Similar  experiments  have  also  been  carried  out  by  Walther 
(1061),  who  in  addition  observes  the  pulse  compression  of  a  frequency  modulated 
acoustic  wave. 

Transient  signals  in  waveguides  have  been  observed  experimentally  by  Saxton 
and  Schmitt  (1963),  and  Ito  (1964,  1965).  The  best  measurements  are  those  by 
Ito  (1965)  who  measures  the  transient  response  of  a  short  Gaussian  envelope 
carrier  pulse.  The  results  show  that  the  theoretical  solutions  of  Forrer  (1958)  do 
not  adequately  predict  the  response.  However,  good  agreement  is  obtained  with  a 
solution  based  on  the  method  of  stationary  phase. 

The  dispersion  of  pulsed  dc  and  rf  signals  in  plasmas  have  been  observed 
by  Schmitt  (1964,  1965).  The  response  of  short  unidirectional  pulses  is  oscillatory 
in  nature  with  a  "ringing"  frequency  that  is  characteristic  of  the  plasma 
frequency.  Such  electromagnetic  pulses  can  therefore  be  used  as  a  diagnostic 
tool. 

2.  WAVE  PACKETS  AND  GROUP  VELOCITY  IN  ISOTROPIC  PLASMAS 

The  simplest  ideas  concerning  the  propagation  of  signals  in  a  dispersive 
medium  involve  the  concept  of  a  wave  packet.  A  wave  packet  is  a  signal  which 
contains  a  narrow  band  of  frequencies  centered  about  a  carrier  frequency.  The 
frequency  spectrum  E(u)  is  then  some  type  of  a  peaked  function  and  is  related  to 
the  time  function  £(t)  by  the  Fourier  transform  pair 


+• 

*(t)  »  /  E(u*  eiu,t  du 

-00 


(1) 


E(u* 


+00 


/ 


F(t)  e‘iwt  dt. 


(2) 

1 
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For  example,  for  a  Gaussian  envelope  carrier  pulse 

_jti 

,~.2  iui  t 

^(t)  -  Eq  e  21  e  0 
the  Fourier  transform  is 


E  +°°  ,T2  t(u>  -  w)t 

EM-»?  /  *  2T  «  ° 


dt 


F  +0°  2 

SW-lf  /  eat‘bt  dt. 

-OO 


where 


a  1  -i  (u  -  uQ) 


h  * 


2T 


2  ‘ 


(3) 


(4) 


Evaluating  the  integral  in  Eq.  (4)  (see  Appendix  A)  and  substituting  from  Eq.  (5) 
one  obtains 


E 

E(w)  «  -2. 

V2I 

This  frequency  spectrum  together  with  the  corresponding  time  response  is  shown 
in  Figure  la. 

Similarly,  for  the  rectangular  envelope  carrier  pulse  shown  in  Figure  lb  the 
frequency  spectrum  is  given  by 


T  exp 


!■ 


T2  (to)  -  u»ri)2 


(6) 


E(w)  «  E0 


sin  (h)  -  a>0)T 

»  (to)  •  to)Q)  ' 


(V) 


Note  that  in  each  of  these  cases  the  frequency  spectrum  becomes  narrower  and 
more  peaked  as  the  width  of  the  time  response  increases.  Thus  the  concept  of  a 
wave  packet  can  be  expected  to  break  down  wher  the  carrier  pulse  becomes  very 
short. 


FIG  lb  RECTANGULAR  ENVELOPE  CARRIER  PULSE 


Figure  1.  Frequency  Spectrums  for  Two  Types  of  Wave  Packets 
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Now  consider  the  propagation  of  a  wave  packet  in  a  dispersive  medium. 
Suppose  that  at  z  *  0  the  time  response  is  given  by^(t,  0)  wiih  the  corresponding 
frequency  spectrum  E(u>,  0).  Then,  if  the  medium  is  lossless,  at  some  point,z, 
each  frequency  component  will  have  undergone  an  appropriate  phase  change  and  the 
resultant  time  response  will  be  given  by 

+00 

^(t,  z)  =  /  E(w,  0)  e1^*1^'  z)  du).  (8) 

-  OO 


Fci  an  isotropic,  lossless,  homogeneous  plasma 

0(o>,  z)  =  k  (w)  z  (9) 

where 


is  the  wave  number  and  n  is  the  electron  plasma  frequency. 

If  E(u>,  0)  is  sufficiently  peaked  about  wQ,  then  0(u>)  can  be  expanded  in  a  Taylor 
scries  about  uQ  and  gives 

0(w)  *  0  (wQ)  +  (w  -  wQ)  $HwQ)  +  £  (w  -  «0)2  0  "(w0)  +.  .  .  ,  (11) 

2 

or,  keeping  only  terms  up  to  order  (<*>  -  wQ) 

♦  *  *0  +  V  +  2  y2  <12> 
where 

V  ■  w  -  w0.  (13) 

Substitution  of  Eq.  (12)  Into  Eq.  (8)  gives 


ea,  z) 


.JKX-  V 


+  00 

/ 


-oo 


E(o>,  0) 


e  i(t  -  *')  V 


(14) 


7 


The  case  of  a  Gaussian  envelope  and  a  rectangular  envelope  carrier  pulse  will  be 
considered  separately. 

2.1  Propagation  of  u  Gaussian  Kn\elope  ('.airier  Pulse 

Substitution  of  the  frequency  spectrum  for  a  Gaussian  envelope  carrier  pulse 
given  by  Eq.  (6)  into  Eq.  (14)  gives 


...  i  (w  t  -  a  ) 

tf(t.  z>  -  Eq  —L  e  0  ° 

°  V2rr 


I  exlJ  .  i(t-P')  V  -  (-y  »  t-f)  V2' 


T  i(wot  ■  <*o)  r°°  1  2  I 

£(t,  z)  =  E  — -  e  J  exp  ivV  -  (IV  d  V' 

°  -/TV  L  I  \ 


^  .  _  T  i(W  /r  ft 

^(t'  Z)  =  Eo  e  e 


(see  Appendix  A),  where 


a  -  i  (t  -  0’) 


0  -  £  <T2  +  i  0"). 


From  Eq.  (16)  it  is  found  that 


4?  *  *>2  +  iC 


where 


y2  .  - - 

-2  >♦«)’ 


>”  (t  -  _ 
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Thus,  the  solution  Eq.  (15)  can  be  written  in  the  form 


(19) 


where 


(20) 


From  Eq.  (19)  it  is  seen  that  the  wave  packet  has  become  frequency  modulated 
while  the  envelope  has  remained  Gaussian  but  has  decreased  in  amplitude  and 
spread  out.  The  maximum  of  the  pulse  envelope  occurs  at  the  delayed  time 


Ik 

d  u) 


—  where 

Vg 


dj£ 

8k 


oo 


(21) 


is  the  group  velocity  of  the  wave  packet.  From  Eq.  (10)  and  Eq.  (21)  the  group 
velocity  for  an  isotropic  plasma  can  be  written  as 


A  plot  of  vjc  vs  n  /w  is  shown  in  Figure  2.  For  this  same  plasma  we  can  write 
f  (w  )  as 

v  ^ 


(23) 


g 


HI 


Figure  2.  The  Group  Velocity  as  a  Function  of  the  Plasma  Frequency 


Figure  3  3hows  the  relation  between  the  normalized  plasma  frequency 


P 


(24) 


and  the  damping  factor  y  |  (wQ)  j  plotted  from  Eq.  (23). 

2.1.1  APPLICATION  TO  PULSE  DEFINITION  AND  SEPARATION  PROBLEM 

The  results  of  the  previous  section  will  now  be  applied  to  the  problem  of  re¬ 
solving  two  Gaussian  pulses  which  are  close  together.  Consider  a  double  Gaussian 
pulse  as  shown  in  Figure  4a.  After  traveling  through  a  plasma  a  certain  distance 
the  two  pulses  may  become  dispersed  to  the  ciegree  shown  in  Figure  4b.  This 
condition,  at  which  the  amplitude  of  each  pulse  is  equal  to  e  0,5  at  (t  -  0')  = 

r*  ,  is  defined  as  the  condition  of  no  pulse  definition.  From  Eq.  (19)  and  Eq.  (17) 

it 

this  condition  occurs  when 
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Figure  3.  Relation  between  Normalized  Plasma  Frequency  and  the  Damping  Factor 


or  when 


"V  1  + 


Equation  (2S)  is  plotted  in  Figure  5  and  separates  the  region  of  good  pulse  definition 

from  the  region  of  no  pulse  definition.  This  figure  indicates  the  critical  ratio  of 

T 

pulse  separation  to  pulse  width  for  different  values  of  ~p==~  .  Values  of 


|  0"  |  may  be  convenientely  determined  from  Figure  3  for  different  plasma 
frequencies  and  propagation  distances.  Due  to  the  quasi -monochromatic  nature  of 
the  wave  packet,  the  pulse  definition  given  by  Eq.  (25)  is  only  a  rough  approxima¬ 
tion.  A  more  accurate  treatment  can  be  carried  out  using  the  methods  of 
saddleooint  integration. 
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Figure  5.  Regions  of  Good  Pulse  Definition  and  No  Pulse  Definition 


2.2  PwpmllH  of  a  Rectaagalar  Eavelepe  Carrier  Paiae 

Consider  now  a  rectangular  envelope  carrier  pulse  which  has  the  frequency 
spectrum  given  by  Eq.  (7).  Using  this  expression  in  Eq.  (14)  one  obtains 


«*.$.**»  /*  atp  dV 


E  e 


i(w0t-^0) 


+  00 
/ 


exp{i(t-$'  +  T)V-i-£y  V2}dV 


+  00 


/ 


d  V 


(26) 
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Eo  e  °  J"  "  expjiaj  V-b  V2}  d  V  *  exp{  i ag  V-b  V2}d  v] 

2  it  i  J  V  J  V 

L.-00  -oo  J 


S(  t,z)  = 


(26) 

(Cont) 


where 


ax  =  f  +  T 

a2  =  *'  *  T 

■  _  i£ 
b  "  2 


t'  =  t  -  0' 


(27) 


The  integrals  in  Eq.  (26)  can  be  evaluated  from  the  results  of  Appendix  B.  One 
then  obtains 


E  i(w  t-0  ) 

^■(t,z)  =  -jp  e  |  erf 


From  Eq.  (27)  it  is  found  that 


al  .  (1  -  i)  (t1  +  T) 

2/7  -T-  VlTT  • 
-4=  tut.  . 


-  erf  I 


[2/1 


'  T'  ViTT 

Substituting  Eq.  (29)  into  Eq.  (28)  and  uaing  the  Identity 
erf  £^(l-i)zj  ■  <1  -  i)  F(z) 


(28) 


(29) 


(30) 


where 


F(z)  ■  /  exp  (i^r1-)  du  ■  C(z)  +  i  S(z) 


(31) 
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is  the  Fresnel  integral  function,  one  obtains 


E_  i(uL :  -  r  -i 

£(t,z)  =  -f  (1  -  i)  e  °  °  |_F(Aj)  -  F(A2)J 

where 


a  •  H-t-I 


The  envelope  of  the  time  response  Eq.  (32)  is  then  given  by 
|*<t,x)|»  +  Y2 


(32) 


(33) 


(34) 


wnere 

X  *  C(At)  -  C(A2)  , 

Y  »  StAj)  -  S(A2)  .  (35) 

♦i 

Knop  and  Cohn  (1963)  have  plotted  this  envelope  aa  a  function  of  Tj-f  for  different 

o-J  |  a"| 

values  of — .  Ginzburg  (1962)  has  plotted  the  envelope  as  a  function  of 

t'  1  2  T 

... .  for  different  values  of  . 

V»vi  V*I7i 

The  discussion  thus  far  has  dealt  with  the  ideas  of  a  wave  packet.  This  notion 
requires  that  the  frequency  spectrum  be  narrow  and  peaked.  This  is  clearly  not 
the  case  during  the  time  interval  in  which  the  pulse  is  building  up,  For  during  this 
time  the  frequency  spectrum  is,  in  general,  very  wide.  Thus,  the  results  of  the 
previous  sections  are  limited  and  would  not  be  expected  to  hold  during  the  initial 
arrival  and  main  buildup  of  the  pulsed  Bignal.  In  order  to  study  this  facet  of  the 
problem,  a  more  general  treatment  of  signal  propagation  in  plasmas  will  now  be 
discussed. 
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3.  GENERAL  ANALYSIS  OK  TRANSIENT  SIGNAL  PROPAGATION 

The  general  method  of  obtaining  the  time  response  of  a  sigi  propagating  in 
an  isotropic,  lossless  plasma  will  be  considered.  The  analysis  will  be  based  on 
Laplace  transform  techniques.  The  general  integral  to  be  evaluated  may  be 
obtained  from  Eqs.  (8),  (9),  and  (10)  by  letting  s  =  i  w  and  picking  up  a  factor  oi  2  77, 
and  is  written  as 


<?(t,z)  =  /  E(s,  0)  exp 

^s 

This  expression,  which  is  derived  ndependently  in  Appendix  C,  is  just  the  inverse 

|  z  JJ~2  I 

Laplace  transform  of  EL  ,  0)  exp  j-  —  ys  +IT  ^  where  E(s,  0)  is  the  Laplace 
transform  of  the  time  response  at  z  =  0.  The  contour  is  a  straight  line  between 
-  i°°and  oo  +  i<*>  where  <jq  is  to  the  right  of  all  singularities  in  the  complex 
s-plane  (s  =  ff  +  i  u). 

There  are  very  few  initial  time  responses  for  which  the  inverse  transform 
given  by  Eq.  (36)  may  be  found  from  tables.  One  such  time  response  whose 
solution  may  be  written  down  is  a  turn-on  Bessel  function  whose  frequency  is  the 
plasma  frequency.  That  is, 

£<t,  0)  *  JQ  (  nt)  U(t).  (37) 


^st  -  —  *fs 

\  r  * 


s2  +  n2*  ds. 

* 


(36) 


where 


|0  t  <  0 

U(t)  - 

(  1  t  >  0 


The  Laplace  transform  of  Eq.  (37)  is 


E(s,  0) 


(38) 


J 

V 
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Equation  (36)  may  then  be  written 


tf(t.z)  *  Jf 


-1 


exp 


Vs2  +  n2 


\/s2  +  n' 


(39) 


which,  from  tables,  is 


tf(t.z)  = 


/  •>„  ("V^4) 


0  <  t  <  z/c 


(40) 


A  similar  solution  is  obtained  when  the  initial  time  response  is  a  step  H  field.  To 
see  this,  note  from  the  transformed  Maxwell  equation 


8  E  (s.z) 
-"ST1-" 


-  8  H  (s.z)  I 


that 


-  Vs2  +  n2  E  (a,z)  a  ~  M0  8  H  (s.z)  , 


from  which 


E  <f,z) 


V?  +  ns 


H  (s,z)  . 


(41) 


Thus  for  an  Initial  step  H  field 


*<t,  0)  *  U(t), 


•i- 


H  (s,  0) 
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so  that,  from  Eq.  (41) 


E  (s,  0) 


(42) 


Since  this  initial  transform  is  oi  the  same  form  as  Eq.  (38),  the  resulting  time 
response  will  be,  from  Eq.  (40), 


^(t,z)  = 


1  '•  ("/*-?) 


0  <  t  <  - 
c 


(43) 


t  >-r 


Other  forms  of  initial  time  responses  which  may  be  of  more  practical 
importance  do  not  result  in  transforms  with  known  inverses.  Therefore,  in  the 
following  sections  several  difierent  techniques  for  evaluating  the  integral  Eq.  (36) 
will  be  investigated. 

* 

3.1  Convolution  I nteginls 

The  Laplace  transform  of  a  function  f(t)  is 


CO 

F(s)  *  X  [  f(t)j  *  /  f(t)  e'8t  dt 

o 


(44) 


and  the  inverse  Laplace  transform  of  F(s)  is 


f(t)  «  jr’1  [F(s>]  *  jjy  /  F(s)estds. 


(45) 


't'he  convolution  theorem  states  that  if  F(s)  is  the  Laplace  transform  of  f(t)  and  G(s) 
is  the  Laplace  transform  of  g(t),  then 

t 

(F(s)G(s)]  «  /  f(t  -  t')  g(t')  dt1. 
o 


(46) 
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Application  of  this  theorem  to  Eq.  (36)  by  using  data  from  tables  (for  example,  see 
Erd^lyi  et  al,  1954),  which  show  that 


As  an  example  of  the  use  of  Eq.  (48)  let  the  time  response  at  z  *  0  be  a  step- 
carrier  sine  wave  given  by 


#(t,0)  «  EQ  sin  t  U(t), 

[i  wt  ] 

Eo  *  U(t>  •  (4 

where  Im  means  "tee  imaginary  part  of".  Using  Eq.  (49)  in  Eq.  (48)  one  obtains 
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Normalizing  all  quantities  to  the  signal  frequency  u>0,  one  may  write  Eq.  (50)  as 


e  (r,  n )  =  im  eq 


U(7-t»  (51) 


where 


r 


P 


(o2) 


The  form  of  the  solution  Eq.  (51)  is  interesting.  First,  propagation  of  the  signal 
front  or  wave  front  proceeds  at  the  speed  of  light  in  free  space.  The  solution  con¬ 
sists  of  two  ter  ms.  The  first  term  represents  the  propagation  of  a  plane  wave  in 
free  space.  The  second  term  represents  the  dispersive  properties  of  the  plasma. 
The  two  terms  combine  to  give  the  total  or  dispersed  wave  that  is  propagating  in 
the  plasma. 

Equation  (51)  has  been  evaluated  numerically  and  several  typical  results  for  the 
overdense  case  are  presented  in  Figure  6. 

In  order  to  investigate  the  effects  of  a  finite  rise-time  on  the  dispersion  of  a 
carrier  pulse,  let 


t  U(t)  , 


/  t\ 

*(t,0)  *  EQ  I  1  •  e  J  sin 

f  (  -oui  t\  i  a)  t  j 

[l »“•  J  *  J 


<?(t,0)  *  Im  E 


U(t>. 


(53) 


Substituting  Eq.  (53)  into  Eq.  (48)  and  using  Eq.  (52)  one  obtains 

riuj1  (p  -/u2  -  n2) 


*  Im  E 


ei(r-rj)  _  e*(o-i)(r-i))  .  P|?  eir  j 


+  Ptj  e 


l  _(ar"i)T  J 

-(or-i)T  /  _ _ L 


(pyi^v) 


du 


vi2  -  ,2 


U(r-i». 


du 


(54) 
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Figure  6.  Propagation  of  a  Sine  Wave  Electric  Field  in  an  Overdense  Plasma 


Equation  (54)  has  boen  evaluated  numerically  and  the  results  are  shown  in 
Figure  7. 

As  a  t'urtner  example  consider  a  step  E  field  turned  on  at  z  *  0,  that  is, 

<?(t,0)  *  Eq  U(t).  f5j> 

Then,  from  Eq.  (48),  the  time  response  at  z  will  be  given  by 


(56) 
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(a)  (b) 


Figure  7.  Envelope  Response  for  the  Propagation  of  a  Sine  Wave  Electric 
Field  with  a  Finite  Rise  lime 


This  expression  has  been  evaluated  numerically  and  typical  cases  are  shown  in 
Figure  8. 

3.2  Series  Solutions  and  Loamel  Functions 

It  is  possible  in  the  case  of  a  lossless,  isotropic  plasma  to  obtain  an  exact 
solution  for  a  step  carrier  sine  wave  input  electric  field.  The  procedure  hinges  on 
several  known  inverse  Laplace  transforms  and  several  identities  involving  Bessel 
functions.  Let  the  time  response  of  the  electric  field  at  z  =  0  be  given  by 

6"(t,0)  *  E_  sin  «)t  U(t)  .  (57) 

o  o 

The  transform  of  Eq.  (57)  is 


E(s,0)  *  jrl#(t,0)]  -  E0  -y . - . -J 


S  UL 


and  the  time  response  at  z  is  then  given  from  Eq.  (36)  by 

I 


E(s,0)  exp 


4^  *  n2  |  . 

I  c  v  ) 


(58) 


f(t,z)  =  X 


1 


(59) 


Figure  8.  Propagation  of  a  Step  Function  Electric  Field  in  a  Plasma 
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If  the  expression  for  E(s,  0)  given  by  Eq.  (58)  is  used  in  Eq.  (59),  the  inverse 
transform  of  the  resulting  expression  is  not  known.  However,  it  has  been  found 
previously  that  the  inverse  transform  of  an  input  step  Bessel  function  with  argument 
IT  t  is  known  (see  Equation  (40)).  This  suggests  the  possibility  of  synthesizing 
sinwQt  from  a  sum  of  Bessel  functions  with  argument  lit.  To  this  end  the  well 
known  identity  is  used  (see  Abramowitz  and  Stegun,  1964) 


sin  (z  cos  6)  *  2  £  (*l)n  J2n+i  (z)  cos  ((2n+l)0i 


n  =  o 


2n+l 


(60) 


By  letting 


z  *  nt 

3i 


cos  e  *  n 

cosh  a  =  “r 


“o  *  n 


(61) 


where  9  -  ia  for  >11,  one  can  write  Eq.  (60)  as 


sinwot»2  £  <*Dn  An  J2n+1  (Ht) 


n  *  o 


(62) 


where 


cos  (2n+l)  0 


"o  *  11 


i.S 

n  i 

(  cosh  (2n+l)  s  «0  £  n 


(63) 


Now,  from  tables 


2n+l 


(64) 


so  that,  from  Eq.  (62)  one  may  write 


jc  [  sin  wQt  U(t)l  •  y==  f  <‘1>n  An  ( - - -) 

Vs2  +  n2n“°  \s  +  v»z  +  n2  / 


2n+l 


(65) 
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One  therefore  obtains  from  Eqs.  (59)  and  (65) 


€(t,z)  » 


(66) 


Since  the  series  in  Eq.  (66)  is  uniformly  convergent,  one  may  interchange  the 
inverse  and  summation  operations.  If  this  is  done,  and  if  the  transform  pair 
(see  Erdllyi,  et  al,  1954), 


.-1 


exp 


{-I  } 


n 


2n :  1 


fs2  +  n2 


<S+  Vs^+n2 


'± 

•♦f 


2n+l 
2 


2n+l 


U(t-|). 


(67) 


is  used  one  may  write  Eq.  (66)  as 


.  .  2n+l  _ 

■  JEo  fo  a„  (n  ^ 


-=5  ]  U(t  (68) 


Normalizing  Eq.  (68)  to  wQ  by  using  the  relations  in  Eq.  (52)  one  obtains 


#(T.l»>  -  2E0  U(t-8)  l  (-1  f  An  ^ )  J2n+, 


(89) 


where 


cos  (2n+l)  6 


I 


cosh(2n+l)  a 


P^l 
P  £  1 


(70) 
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and 


COS  0  s  p" 

P  >  1. 

cosh  a  *  p" 

P  <  1. 

(71) 

It  is  seen  then  that  the  electric  field  response  for  this  case  can  be  represented  as 
a  convergent  infinite  series  in  terms  of  Bessel  functions.  This  result  has  been 
obtained  by  Knop  (1964).  Similar  results  have  been  obtr.ned  by  Cerrillo  (1948), 
Rubinowicz  (1950),  and  Gajewski  (1955). 

Foi’  the  underdense  case  of  propagation  it  is  possible  to  write  Eq.  (69)  as 
the  sum  of  two  Lommel  functions  of  order  one.  The  Lommel  function  of  the  first 
kind  (w,  z)  is  defined  by  Dekanosidze  (1960)  as 


2n+i> 


U„  (w*  z)  *  l  (  l)n  <*>  J2n+„(z)- 


(72) 


n  =  o 


For  the  underdense  case, P  <,  1,  the  An  in  Eq.  (69)  is  equal  to  cosh  (2n  +  l)a, 
where 


a  =  cosh  1  (^) 


(73) 


or,  noting  that 

ea  »  cosh  a  +  sinh  a, 
one  may  write 


If  one  lets 

a  * 


(74) 


then 


a  *  In  a 
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and 


(2n  +  l)  a  *  In 


J2n+l) 


(75) 


Using  Eq.  (75)  one  can  write 
cosh  (2n  +  1)  a  =  j- 


(2n+l)a  _ i _ 

g(2n+l)a 


or 


cosh  (2n  +  l)a 


y  exp  In 


(2n+l) 


exp 


Jlno<2n+1>}J 


and  finally 


(20+1)  -(2n+l) 

An  1  cosh  (2n  +  l)a  *  - g —  +  3 - g —  . 


The  general  term  of  the  series  in  Eq.  (63)  then  becomes 

(2n+l) 


(76) 


|  frit 


Using  the  identity 

.  p  iTS.nL 

VrM)  /j — ? 

PVT  -  «r 


one  can  write  the  general  term  as 

*Ht»n 


H. 

/  P(r-q)  \an+1  /  _ \1 


(77) 
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One  defines  the  following  quantities 


Wj  =  a  P(r-f)) 


(r-tj) 


w2  •  -g(T-tj)  =  f '  ? 


[if  l/l  -  P2  ] 


(r-t?) 


2  =  P 


n  2 

VT  -  rr 


(78) 


where  Eq.  (74)  has  been  used.  Using  Eqs.  (69),  (72),  (78)  and  Eq.  (77)  for  the 
general  term  of  the  series  one  may  write  the  solution  as 

t  (T,r\)  •■=  Eq  U(r-Tl)  Pui(w1,z)+  (w2,  z)  j  .  (79) 


Lommel  functions  of  two  variables  have  been  tabulated  by  Dekaoosidze  (1960)  so 
that  Eq.  (79)  may  be  ”eadily  evaluated  for  a  given  P,  »)  and  r. 

In  a  similar  way  it  may  be  shown  that  Eq.  (79)  is  also  the  solution  for  the 
overdense  case,  P  >  1,  if  w^  and  w2  are  defined  by  the  new  relations 


«  [i  +  iy/p2  -  1  ]  <T-t )), 


w2  = 


[l  +  iVp2  - 1  ] 


(r-l)). 


(80) 


The  arguments  of  the  Lommel  functions  are  complex  for  this  case  and  since 
these  are  not  tabulated  the  usefulness  of  this  representation  is  somewhat  limited. 


3.3  (aatsw  latsgrallMs 

The  solution  to  the  problem  of  transient  signal  propagation  in  an  isotropic, 
lossless  plasma  is  given  by  Eq.  (36),  namely 


£(t,z) 


J  E(s,0)  exp  |  st  y42  +  P"  J 

y  I  I 

s 


ds. 


(31) 


where  the  contour  yg  lies  to  the  right  of  all  singularities  as  shown  in  Figure  9a.  If 
the  contour  is  closed  in  a  large  semi-circle  of  radius  R  to  the  right  (see  Figure  9b) 
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Figure  9.  The  Complex  Frequency  Plane 


so  that  s  >  >  II,  then  the  exponent  in  Eq.  (81)  becomes  s(t-z/c)  and  for  this  case 
the  integral  along  this  large  semi-circle  goe3  to  zero  as  R  — *  «  when  t  <  z/c. 
Since  the  integral  around  the  closed  contour  in  Figure  9b  is  zero  by  Cauchy's 
integral  theorem,  ther.  the  integral  Eq.  (81)  will  be  zero,  so  that 

<r(t.z)  *  0  i  <|  (82) 


as  has  beat  found  previously.  Thus  no  signal  can  arrive  at  a  point  z  with  a  speed 
g:  eater  than  that  of  light  in  free  apace. 

For  t  >  z/c  the  contour  will  have  to  be  closed  in  a  large  circular  arc  to  the 
left  in  order  to  make  the  contribution  along  this  circle  vanish.  However,  when  this 
is  done  singularities  will  be  enclosed  in  the  complex  plane  which  will  contribute  to 
the  value  of  the  integral  in  Eq.  (81). 

Two  types  of  singularities  are  generally  encountered.  The  first  are  poles  of 
the  function  E(s,  0)  which  will  give  rise  to  the  steady  state  solution.  The  second 
are  branch  points  which  arise  from  the  double-valued  nature  of  the  function 


Vf  o  o 

s  ‘  +  FT  .  The  branch  points  of  this  function  occur  at  tin  and  in  order  to  carry 
out  the  integral  in  Eq.  (81)  the  complex  plane  must  be  cut  in  such  a  way  as  to  make 


2!* 


/  2  2 

the  function  vs  +n  single-valued.  It  will  be  found  that  the  integration  along  the 
branch  cuts  gives  rise  to  the  transient  solution  to  the  propagation  problem.  Seveial 
ways  o:  'aking  the  branch  cuts  are  possible  and  two  ol  these  ways  will  ue  investigated 
in  detail. 

3.3.!  BRANCH  CUT  ALONG  IMAGINARY  AXIS  FROM  -in  TO  ^in 

The  fi r st  step  in  evaluating  the  integral  in  Fq.  (81 )  is  to  study  the  double  valued 
function 


f(s)  - 


(83) 


One  wishes  to  form  two  Riemann  sheets  by  cutting  the  complex  s -plane  along  tin- 
imaginary  axis  from  -in  to  rill.  In  order  to  properly  define  these  two  sheets  one 
considers  f(s)  to  be  the  product  of  the  two  double-valued  functions 


and 


Vi 


s  *  in 


(8-U 


F or  each  of  these  functions  one  takes  the  branch  cuts  shown  in  Figure  10.  Sheet  l 
of  Vs  -  ifl  is  defined  in  the  interval  ^  and  Sheet  2  in  the  interval 

*£  <  $i<  •  Sheet  l  of  Vu  +  in  is  defined  in  the  interval  -y  <  o.,  <  , 


3  1  V 

w'hile  Sheet  2  is  defined  in  the  interval  -y  <  ^  <  —y  .  Using  these  branch  cuts, 

one  k<oks  at  the  product  of  the  functions  V »  -  if!  end  V*  +  to  determine  the 

/2  2 

ign  distribution  for  the  function  f(s)  *  Vs  +  n  .  The  Riemann  sheets  for  f(s) 
must  be  defined  with  respect  to  both  ^  and  ^  Sheet  1  of  f(s)  is  defined  by  the 
intervals  of  ^  and  02  given  by 


1 


< 


3  P 
2 


3  V 
2 
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Figure  10.  The  Complex  Frequency  Plane  for  a  Branch  Cut  from  —in  to  +  ifl 


which  are  obtained  from  the  product  of  Sheet  1  of  v4*  -  in  and  Sheet  1  of  V  s  -*  in  . 
These  conditions  are  identical  to  the  intervals 


< 


7  I 

2 


if 

2 


<*2  < 


7  * 
2 


which  are  obtained  from  the  product  of  Sheet  2  of  Vs  -  in  and  Sheet  2  of  Vs  *  in . 
The  second  shoet  or  Sheet  2  of  f(s)  is  given  by  the  intervals 


n 

2 


< 


3  n 
2 


< 


ft. 


in 


or  by  the  intervals 
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SHEET  I  SHEET  2 

Figure  11.  The  Sign  Distribution  on  *.ach  Riemann  Sheet  for  the  Branch  Cut 
of  Figure  10 
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The  integration  on  Sheet  1  will  now  be  carried  out  for  a  step  sine  wave  input 
along  the  contour  shown  in  Figure  12.  For  this  case 


E(s,  0)  E  -2  9-  ,  (86) 

8  +Wo 


Figure  12.  The  Contour  for  the  Branch  Cut  in  Figure  10 


One  may  readily  show  that  as  R  —■*  *>  the  integral  along  the  large  circular  arc  C 
vanishes.  The  two  integrations  along  the  negative  real  axis  cancel  and  as  the 
radius  of  the  two  small  circles  c+  and  c#  around  each  branch  point  goes  to  zero 
the  contribution  to  the  integral  from  these  small  circles  also  vanishes.  Thus  the 
only  contributions  to  the  integral  in  Eq.  (81)  come  from  the  residues  of  the  poles  and 
from  the  integration  along  each  side  of  the  branch  cut.  Therefore 

f  +  J  *  2  si  2  Residues.  (87) 

yB  Branch 
Cuts 
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Along  the  branch  cut 


Px  =  n-  w 
p2  =  n  +  w  . 


(88) 


Using  Eq.  (88)  in  Eq.  (85)  together  with  the  sign  distribution  of  Sheet  1  in  Figure  11 

/  2  2 

one  sees  that  on  the  right  side  of  the  branch  cut  f(s)  *  +  V  fl  -  cu  while  on  the 


left  side  f(s)  * 


-  J fl2  -  u>2  . 


From  Eqs.  (81),  (86)  and  (87),  one  can  then  write 


...  ,  Eowo  riTI  exp  {  st  |  t/n2  -  tu2  } 

£U,Z'  2  w  i  J 

-in 


2  2 
s  +  cu 
o 


ds 


r111  exp  {  st  -  ~  Jn 2  -  u>2  }  v. 

__2.  j  - S - - - ds  =  2  ffi  Residues, 

+in 


2  2 
3  +  <*> 
o 


or,  letting  s  =  i  cu 


-  e 


Ecu  +n  eiu»t  f  |Vn2  -  «2  -U 

...  ,.oor  e  I.  c  »  c» 

f (t.z)  +  -j-j-  J  ;-r  5  e 

-n  wo  "w  \ 

3  2fi  T]  Residues, 


du 


which  becomes,  after  evaluating  the  residues. 


g(t.z)  *  Eq  sin 


E  w  eiwt8inhfv/n2-w2 


I  *  ta,°z  faik  rn  eiwt,,in 


du). 

(89) 


The  integrand  of  the  imaginary  part  of  the  last  term  in  Eq.  (89)  is  an  odd  function 
and  therefore  integrates  to  zero.  One  can  then  write  Eq.  (89)  in  normalized  form, 
using  Eq.  (52),  as 


e(T,Ti)  *  EQ  sin 


P  f~ 5  2 

j  cos  u  t  ainh  ^  V  P  ■  u  d|] 


-P 


1  -  u 


(90) 
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It  is  evident  that  the  total  solution  is  the  sum  of  a  steady  state  term  and  a 
transient  term.  The  integral  representing  the  transient  term  is  not  particularly 
suited  for  numerical  integration  since  the  integrand  oscillates  rapidly  for  large  T. 
In  order  to  overcome  this  difficulty  the  solution  based  on  a  different  branch  cut 
will  now  be  examined. 

3.3.2  BRANCH  CUTS  IN  LEFT  HALF  PLANE 

A  different  branch  cut  which  better  lends  itself  to  numerical  integration  will 
now  be  considered.  If  the  s-plane  is  cut  as  shown  in  Figure  13  then  the  integra¬ 
tions  along  tbe  branch  cuts  will  be  over  negative  values  of  the  real  part  of  s.  One 
might  expect  the  resulting  integrands  not  to  oscillate  as  much  as  if  one  integrated 
along  the  imaginary  axis. 


Figure  13.  The  Contour  for  Branch  Cuts  in  the  Left  Half  Frequency  Plane 

Let  the  initial  time  response  be  an  exponential  function  of  the  form  E(Jei  wol 
and  take  the  imaginary  part  of  the  result  to  give  the  response  to  a  sine  wave  input. 
The  integral  to  be  evaluated  may  then  be  written  from  Eq.  (81)  as 
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f  exp  {  st  -  |  -J&2  +  n2  } 

I - irrfs - 1 

o 


The  integration  will  be  carried  out  on  Sheet  1  of  Figure  13,  defined  by  the 
relations 


-  in  =  px  e 


+  in  =  p,  e 


-  ”  <  0!  <  * 


-  V  <  <?,  <  It 


One  must  now  determine  the  value  of 


/ 2  ,  _2 
(s  +  n 


s  -  in  A/s  +  in  =  -v/ P,  #>o  e 


2  (*1  +  ^2) 


for  each  of  the  four  sides  of  the  branch  cuts  in  Figure  13. 

Consider  first  Side  ©  of  the  upper  branch  cut.  Along  this  side  f  ^  =  +  jr  and 

02  =  "f  +  where  \)/  *  tan*1  m  .  One  then  obtains  from  Eq.  (93) 

•'**•'* 

•  (*1  +  i)  (  co*"^+ 

[-(cos  +  l(co.4-  sin|)] 


+n2 


(-  a  +  i  b) 


[Side  ©  ] 


where 


a  -  cos  +  sin 
b  a  cos  -i^  -  sin 
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Since  0  ^  ^  <  y ,  so  that  0  £  ^  it  is  seen  that  a  >  0  and  b  >  0. 

In  a  similar  manner  one  can  readily  show  that  for  Sides  ©  ,  ©  and  © 
one  obtains 


t/s2  +  n2 

‘V-7-2  <•-*»>) 

[Side 

© 

] 

./2  +  n2 

a  ~  (a  4-  i  b) 

[Side 

© 

] 

-\42  +  n2 

/p 1  P2 

(-»  -  i  b) 

[  Side 

© 

i- 

(96) 

For  each  of  the  four  sides  it  is  also  noted  that 

IoI'n/IjI2  +  4ii  2 
2 


(97) 


Consider  again  Side  ©  .  Along  this  side  s  *  o+ill  and  the  integral  in  Eq.  (91) 

becomes,  using  Eq.  (94) 


Tin?  T+Tn-lu  exp  {<ff+  *nn (“a  +  i b)} 


do 


int 


2»T“  /  -o+  i  (n  -  w  )  e3£P  -  at|da 


(98) 


where 


A  ■  a 


B  »  b 


(99) 


Similarly,  one  con  evaluate  the  integral  along  the  other  three  sides  of  the  branch 
cuts  and  obtain  for  the  total  contribution  from  the  branch  cuts 
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J  *  /  *  / 

Branch  ®  © 

Cuts 


I  ♦  I 
CD  CD 


Noting  that 


_ 1 _ 

'  0  +  i  (11  +  «o) 


+  a-  i  (n+  w  ) 
a2  +  (II  +  wc)2 


(101) 


(102) 


(103) 
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and 


a  *  cos  B  3inh  A 

P  *  sin  B  cosh  A  .  (104) 


As  in  the  previous  section  it  can  readily  be  shown  that  the  contributions  to  the 
integral  from  the  large  circular  arc  in  left -half  plane  and  from  the  small  circles 
around  the  branch  points  vanish,  so  that  the  total  response  to  a  sine  wave  input  is 
found  from  the  relation 

f(t,z)  +  Im  /  =  Im 

Branch 
Cuts 


£2  wi  £  Residues 


or 


\  r  /  1 - v  ■ 

1 

)eo  exp  I  i  (Wot  yw02 -n2 j 

-  / 

1  L  '  '  J 

Branch  1 
Cuts  , 

(105) 


Normalizing  all  quantities  to  <t)Q  according  to  Eq.  (52)  one  obtains  by  substituting 
Eq.  (102)  into  Eq.  (105) 


#(r,ij)  *  Eq  {sin 


where 

M  *  t*  ai  "  V 
N  -J-  (I3  -  1,1 

0  •  tan"1 


N2  sin  £pt  +  0  (r)J  j. 


(106) 


(107) 


and  where  the  integrals  Ij,  Ij,  I3,  and  are  given  by  the  following  relations 


I 
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I.  -  /  aY^P+  V  e'TX  dx 
‘  o  x2  +  (p+ir 


I  =  /  .SUL±§{P-l)  e*TX  dx 
^  o  x^  +  (p  -  ir 


I,  =  f  A*zJLiPt  1)  e"TX  dx 
J  o  x^  +  (p+ir 


*4  =  / 


/3x  -  a  (P  -  1)  _-tx 
— 2 — ■ — ''■■■ . —"i»  •  e  dx 


+  (P  -  1)A 


where 


or  = 

fi  • 

A  = 

B  * 
*  - 


cos  B  sinh  A 
sin  B  cosh  A 


V  |cos-^  +  sin  j 

T)  |cos^  -  sin  ^  j 
tan"1  (j^,)  • 


(108) 


(109) 


The  integrands  in  Eq.  (108)  decay  exponentially  for  large  t  and  for  large  x. 
They  therefore  lend  themselves  to  numerical  integration  since  the  upper  limit  of 
integration  can  be  cut  off  at  some  large  value  of  x  without  appreciable  error.  These 
integrals  were  evaluated  numerically  by  cutting  off  the  upper  limit  at  a  value 
x  *  15/(  t  ■  1).  The  results  are  shown  in  Figures  14,  15,  and  16.  Figure  14  is  a  plot 

of  the  transient  envelope  v42  ♦  N2  as  a  function  of  r.  The  total  transient  solution 

■\A*a  +  N2  sin  [  Pr  +  0]  is  plotted  in  Figure  15.  Note  that  it  is  a  phase  modulated 
oscillation  about  the  plasma  frequency.  The  total  response  found  by  adding  Figure 
15  to  the  steady  state  solution  is  shown  in  Figure  16. 

All  of  the  numerical  results  discussed  thus  far  in  this  report  lose  their  useful¬ 
ness  for  very  large  values  of  I ).  This  is  due  to  the  fact  that  under  these  conditions 
it  is  necessary  to  take  the  difference  of  two  very  large  numbers  in  order  to  obtain 
a  small  number.  Thus  the  errors  can  become  very  large.  However,  it  is  just  in 
this  region  of  very  large  t)  that  asymptotic  solutions  to  the  transient  propagation 
ar^blem  work  very  well.  Asymptotic  solutions  to  the  transient  wave  propagation 
problem  will  be  considered  in  a  second  and  subsequent  report. 


40 


(a) 


4.  CONCLUSION 

i  2  the  propagation  of  wave  packets  In  dispersive  media  was  discussed. 

n.  lilted  i.  .  ».<»>  co.c.p.  “  “-1"  ™  “  deK“  ‘T.'Te 

Th.  grotip^v.locity  l.nnl.ion  l.  «l»«>»-  Eq.  Oil.  T.0  W«  phOh.-.^ 

G.uAlBP^env.lop.  carrier  »d  Ih.  .,u...  P-**  crier.  ««r.  inv... lg.«d.  Th« 

„.v.  p.ch.l.  .or.  loon.  <0  prop.*™  m  , 

second  derivative  of  the  wave  number  with  respect  to  frequency  g  e 
distortion  of  the  wave  packet.  This  distortion  places  an  upper  limit  on i  the _  maxt- 
m  rate  at  which  repetitive  pulses  can  be  transmitted  and  distinguished  in 
mUm  rBtC  ^  "1  ATrude  method  for  estimating  this  maximum  rate  for  Gaussian 

...  P~  -  — >  - — 2  - — 

to  quasi -monochromatic  signals. 


I 
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(«) 


Figure  15.  The  Total  Tranr*  unt  Solution  Obtained  by  Integration  Along  the 
Branch  Cuts 


In  Section  3  a  more  exact  treatment  of  the  propagation  of  transient  waves  in 
isotropic  plasmas  wa«  given.  Laplace  transform  methods  were  used,  and  the 
problem  waB  reduced  to  evaluating  an  inverse  Laplace  transform.  The  problem 
of  a  sine  wave  electric  i.  turned  on  at  t  »  0  in  an  isotropic  plasma  was  investi¬ 
gated  in  detail.  The  solution  can  be  expressed  as  a  series  of  Bessel  functions,  and 
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to 


Figure  16.  The  Total  Solution  of  a  Propagating  Sine  Wave  Electric  Field 


this  solution  can  be  expressed  as  the  sum  of  two  Lommtl  functions.  The  solution 
cast  also  be  given  as  a  convolution  integral  representation.  By  using  methods  oi 
contour  integration,  still  other  integral  representations  can  be  obtained,  and  two 
such  representations  were  derived.  Results  from  numerical  evaluation  of  the 
Integral  solutions  were  presented.  Due  to  limitations  in  the  numerical  methods  cf 
integration,  the  integral  solutions  are  useful  only  for  relatively  small  values  of 
the  parameter  rj .  For  larger  values  of  rj,  asymptotic  solutions  can  be  obtained, 
and  this  will  be  the  subject  of  a  subsequent  report. 
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Appendix  A 


Evaluation  el  I 


1 


u> 


at  -  kt 


dt 


Consider  the  integral 


/ 


at-bt 


dt. 


(Al) 


If  one  completes  the  square  in  the  exponent  of  Eq.  (Al)  and  makes  the  change  of 
variables  y  *  vz  (t  -  -jj~)  one  obtains 


from  which 


.  n 

h  'VT  exP  1  -4b  • 


(A2) 


B1 


Appendix  B 


Evaluation 


1  f  t*p|ioV-bV2) 

of  i, - J  — ! - 'd 

2  J  v 


Consider  the  integral 


l2  *  f 


iiaV  -  b  V2} 


Letting  V  «  ^  and  X  *  one  can  rewrite  Eq.  (Bl)  as 
r  .  -1-  f  exp  jiXu-u2}  . 


l2  *  2il  / 


Noting  that 


f  «lUX  -«  “  ^ 


B2 


so  that 


eiuX  .  1 


—  +  i  / 

u  u  J  _ 


iux 

e  dx 


one  can  write  Eq.  (132)  in  the  iorm 


'  -  — L-  f 
‘2  2  Vi  J 


00  -u2 

r  e 


X  *> 

du  +  -jy  J  dx  /  exp  j  ixu  -  u“  J  du,  (F33) 


where  the  order  of  integration  has  been  interchanged  in  the  second  term. 

The  first  integral  in  Eq.  (B3)  can  be  evaluated  by  using  Cauchy's  integral 
theorem  and  integrating  around  the  pole  at  u  =  0.  One  readily  finds  the  value  of 
the  integral  to  be  equal  to  1/2.  In  the  second  integral  in  Eq.  (B3)  the  integral  over 

du  is  equal  to  VF  exp  ^--2—-^  frorn  the  results  of  Appendix  A.  Therefore, 


l2  =  2  +  2VW 


Making  the  change  of  variable  w  *  x/2  and  using  the  definition  of  the  error  function 


•(f)-  * 


•w 

e  dw 


one  obtains 


I2  "  +  ♦ 

from  which,  since  X  *  a/VTT, 


I2  *  7  1  +  erf 


(l/f)] 


Cl 


Appendix  C 

Derivation  of  Integra!  Expression  (or  tf(t,i) 


Maxwell's  equations  which  describe  the  propagation  of  electromagnetic  waves 
in  an  isotropic  plasma  are 

curl  e  -  -  ^  77- 

96 

curl  H  *  .£  4  fo  "af  (C1) 

where 

J  *  -  Ne  v  .  (C2) 


N  is  the  electron  number  density  and  v  is  determined  from  the  equation  of  motio» 


(C3) 


Consider  the  one-dimensional  problem  in  which  £  is  linearly  polarized  in  the 
x-direction  and  is  propagating  in  the  z-direction.  The  x-component  of  £  (t,  z)  will 
be  written  as  z)  and  the  Laplace  transform  of  6  (t,  z)  as  E(s,  z).  If  one  then 
takes  a  Laplace  transform  in  time  of  Eqs.  (Cl),  (C2)  and  (C3)  and  solves  for 
E(s,  z)  setting  all  initial  conditions  equal  to  zero,  one  readily  obtains  the  equation 


C2 


a2  E(s,  z)  i 
o  z  c 


(s2  +  n2)  E(s,  z) 


0. 


(C4) 


The  object  is  to  determine  the  time  response  «(t,  z)  in  the  semi-infinite  region 
z  >  0  when  the  time  response  e(t,  0)  is  prescribed  at  z  =  C.  The  solution  of 
Eq.  (C4)  may  then  be  written  as 


E(s,  z)  =  E(s,  0)  exp 


(C'5) 


where  E(s,  0)  is  the  Laplace  transform  of  f(t,  0).  The  time  response  z)  is  then 
obtained  by  taking  the  inverse  of  Eq.  (C5).  That  is. 


ff(t,  z)  =  ji-r  /  E(s,  0)  exp  ■  st  - 


s2  +  n2 


ds. 


(C6) 


which  is  Eq.  (36)  in  Section  3  of  the  text. 
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